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ABSTRACT 

We present a precise photometric calibration of the first 1.5 years of science imaging from the Pan- 
STARRS1 survey (PS1), an ongoing optical survey of the entire sky north of declination —30° in 
five bands. Building on the techniques employed by Padmanabhan "et all J2008) in the Sloan Digital 
Sky Survey (SDSS), we use repeat PS1 observations of stars to perform the relative calibration of 
PS1 in each of its five bands, solving simultaneously for the system throughput, the atmospheric 
transparency, and the large-scale detector flat field. Both internal consistency tests and comparison 
against the SDSS indicate that we achieve relative precision of < 10 mmag in g, r, and ipi, and ~ 10 
mmag in z and yp\. The spatial structure of the differences with the SDSS indicates that errors in 
both the PS1 and SDSS photometric calibration contribute similarly to the differences. The analysis 
suggests that both the PS1 system and the Haleakala site will enable < 1% photometry over much of 
the sky. 

Subject headings: Surveys: Pan-STARRSl 
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1. INTRODUCTION 

A central problem in astronomy is relating the number 
of photons recorded at a detector to the physical flux of 
photons emitted from a source. This relation depends 
on important astrophysical parameters, like the distance 
and Galactic extinction to the source, as well as more 
ephemeral, local phenomena like the weather at the tele- 
scope and the sensitivity of the detector. The problem 
of photometric calibration is to characterize these latter 
phenomena to render more universal the detected astro- 
nomical phenomena. 

The photometric calibration of optical data is often 
performed by comparing multiple observations of the 
same sources for large sets of sources and demand- 
ing that their fluxes not change over time. This is 
the same technique used to calibrate cosmic microwave 
background and radio data, and it ha s been exten- 
sively u sed in opti c al ast ro nomy, e.g. [Maddox et all 
(119901). iHonemj t tl (fl992h . iFong et all (Il992l jl99l . 
IGlazebrook et all (|1994t ). and Magnie r et all (119921) . In 
the Sloan Digital Sk y Surv ey lYork et all (|2000l SDSS). 
iPadmanabhan et all ([20081 ) applied this technique and 
achieved a photometric calibration accurate at the 1% 
level. 

Often surveys and observations have been calibrated 
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using repeat observations of a small number of fields 
of standard stars, e.g ., the 2 Micron All Sky Survey 
([Skrutskie et all 12006). These calibration observations 
determine the relation between flux and photon-count 
when the fields are observed, and the relation during 
science observations is extrapolated from them. In the 
SDSS, however, each set of observations was made to 
slightly overlap other observations, and the network of 
these overlaps was used to simultaneously calibrate all 
of the SDSS observations. Upcoming optical wide-area 
surveys like PS1, DES, and LSST plan to improve upon 
this technique by tremendously increasing the number 
of multiply-observed stars; each plans to image their en- 
tire survey area several times. The dense overlapping 
regions of these upcoming surveys should yield a much 
more tightly constrained photometric calibration. 

Current surveys require photometric calibration as ac- 
curate as possible, ideally to better than the percent 
level. Typical photometric uncertainties from point- 
spread- function modeling reach the 1% level, and absent 
an equally good photometric calibration, calibration er- 
rors will dominate this uncertainty. Additionally, the 
width of the s tellar locus i s abo ut 1% in certain color 
combinations ([Ivezic "eFaTI [2007h . requiring an equally 
good calibration to allow stars to be photometrically 
identified most accurately. Studies of the interstellar dust 
at high Galactic latitudes can be still more demanding; 
photometri c calibration dominates the unc ertainty in the 
analysis of Schla fly and Finkbeinerl ([20111 ). Uncertainty 
in the photometric calibration even at the 1% level can 
contribute to significant variation in the number densities 
of gala xies used for cl ustering studies on large angular 
scales ([Ross et all 120121 ). In short, a number of current 
science projects are limited by photometric calibration 
accuracy. 

In this work, we desc ribe our application of the 
IPadm anabh an et all ([20081 ) photometric calibration al- 
gorithm to the first 1.5 years of PS1 survey data to 
achieve a 1% calibration. This paper is organized as fol- 
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TABLE 1 

Pan-STARRSI Bandbass Parameters 



Filter 


A e ff 


A B 


Ar 


ZP 


fi 


9P1 


481 


414 


551 


24.41 


21.92 


rpi 


617 


550 


689 


24.68 


20.83 


ipi 


752 


690 


819 


24.56 


19.79 


Zpi 


866 


818 


922 


24.22 


19.24 


ypi 


962 


918 


1001 


23.24 


18.24 



Note. — Pan-STARRSI bandpass parameters. The column 
A e ff gives the effective wavelength of each filter in nm, while the 
columns Ab and Ar give the filter blue and red cutoffs in nm. 
The AB zero points are given by ZP, and observed sky bright- 
ness in ma gnitudes per square arcsecond by /j. All values are 
from FTbnrv et a l. (2012), except for the zero points. Thes e are 
marginally discrepant from the values of Tonrv et al. (2012) due 
to the variation in throughput over the PS1 focal plane, and are 
only intended to serve as a rough guide. 

lows: in Section 31. 11 we describe the PS1 survey and 
its current status. Section $2] describes the photometric 
calibration algorithm and its application to PS1 data. 
Section Ogives the results of our calibration of the PS1 
data, and the results of the tests used to verify the cali- 
bration. In §HJ we discuss the stability of the PS1 system 
and the atmosphere in light of these results. Finally, in 
Section ^5] we summarize, mention prospects for the fu- 
ture, and conclude. 

1.1. The Pan-STARRSI System and Surveys 

The Pan-STARRSI system is a w ide-field optical im - 
ager devoted to survey operations (Kais e~et a/.l l2010l ). 
The telescope has a 1.8 meter diameter p rimary mirror, 
locat ed on the peak of Haleakala on Maui (jHodapp et all 
2004). The site and optics deliver a point spread function 
with a full-width at half-maximum (FWHM) of about 
one arcsecond, over a seven square degree field of view. 
The focal plane of the telescope is equipped with the Gi- 
gapixel Camera 1 (GPC1), an array of 60 4800 x 4800 or- 
thogonal transfer array (OTA) CCDs (jTonrv and Onakal 
2009: lOnaka 7TZll2008h . Each OTA CCD is further sub- 
divided into an 8 x 8 array of independently-addressable 
detector regions which are individually read out by the 
camera electronics through their own on-chip amplifier. 

Most of the PS1 observing time is dedicated to two 
surveys: the 3tt survey, a survey of the entire sky north 
of declination —30°, and the medium-deep (MD) survey, 
a deepe r, many-epoch survey of 10 fields, each 7 deg 2 
in size (Cha mbersllin preparation! ). Each survey is con- 
ducted in five broadband filters, denoted gpi, rpi, zpi, 
zpi, and 2/pi, that together span 400-1000 nm. These 
filters are similar to those used in the SDSS, except the 
gpi filter extends 20 nm redward of #sdss while the zpi 
filter is cut off at 920 nm. The ypi filter covers the re- 
gion from 920nm to 1030nm with the red limit largely 
determined by the transparency of the silicon in the de- 
tector. These filters and their absolute cali bration in the 
context of Pan-STARRSI are described in iStubbs et all 
(|2010D and lTonrv eTal\ (|2012l ). The filter bandpasses are 
summarized in Table [TJ 

The PS1 images are processed by the Pan-STA RRSI 
Image Processing Pipeline (IPP) (Magnier 2006). This 
pipeline performs a utomatic bias sub traction, flat 
fielding, astrom etry (Ma gnier et all [2008). photometry 
(Magnier 2007|), and image stacking and differencing for 



every image taken by the system. The approximately 
one trillion pixels per night are processed in a massively 
parallel fashion at the Maui High Performance Computer 
Center. 

The 3n survey is executed so that each time a patch of 
sky is visited, it is observed for about 40 seconds twice, 
at times separated by an interval of about 15 minutes 
(jChambersllin preparation! ). The two observations make 
a transit-time- interval (TTI) pair. These observations 
are used primarily to search for high proper-motion So- 
lar System objects. Each year, the field is then observed 
a second time in the same filter with an additional TTI 
pair of images, making for four images of each part of 
the sky, per year, in each of the 5 PS1 filters. The MD 
observations consist of 8, much longer, ~ 200 second ex- 
posures, dithered in both position and position angle. 

The data used in this analysis were taken primarily 
between February 12, 2010 and June 19, 2011, though 
a small amount of data from as early as June 20, 2009 
is included. Figure [T] shows the number of times each 
part of the sky was observed during this period. The 
left panel gives the total number of times the sky was 
observed, while the right panel gives the number of times 
the sky was observed in photometric conditions ( §2.3j) . 
The median number of times each part of the sky was 
observed in photometric weather is 4 in each band. This 
makes for two independent TTI pairs of observations, on 
average, of each part of the sky, though Figure [T] makes 
clear that this coverage is variable, and the sky around 
right ascension 100° is covered only by one TTI pair of 
observations or not at all in gripi. The MD fields have 
been observed much more often; 100-300 times in grizpi, 
and about 100 times in ypi. 

2. METHODS 

Like all CCD-equipped telescopes, PS1 ultimately 
records the number of photons received from objects it 
targetfl The number of photons recorded depends on: 

1. the transparency of the night sky toward the ob- 
ject, 

2. the throughput of the detector, filter, and optics, 
and 

3. the size and reflectivity of the telescope mirror. 

The object of the photometric calibration is to convert 
the measurements of the number of photons recorded by 
the system to measurements of the incident flux from 
the object, eliminating the signatures of the instrument 
and atmosphere. The calibration ca n be divided into two 
separa te procedures, as described in lPadmanabhan et ali 
(2008): the relative photometric calibration, in which 
the differences in system throughput from observation 
to observation are removed, and the absolute photomet- 
ric calibration, in which the number of photons recorded 
for some particular configuration of the telescope is con- 
verted to a magnitude on the AB magnitud e system, 
which is based on physical units of flux (jOke and Gunnl 

9 The gain of the GPC1 camera is nearly 1 ADU/electron, and 
we include the quantum efficiency (electrons per photon) in the 
detector throughput, so we do not distinguish between ADUs and 
photons in this discussion 
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Fig. 1. — Map o f the number of times the sky has been observed by Pan-STARRSl, overall (left panels) and in photometric weather 
(right panels) (see flOT) . The x-axes give right ascension and the y-axes give declination, both in degrees. White to black spans to 16. 
The ten black circles are the locations of the 10 medium-deep fields, which are observed more frequently than the rest of the survey area. 
In the winter of 2010, the combination of bad weather and the malfunction of the Pan-STARRSl shutter suspended operations, leading to 
an area of poor coverage in gri around right ascension 100° . 
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|1983[ ). This paper presents the relative calibration of 
the Pan- ST ARRS 1 system; th e absolute calibration is 
described in lTonrv et all (j2012[ ). Given an absolute flux 
calibration for a single star, the relative calibration trans- 
fers this absolute calibration over every observation of the 
survey. 

In this section, we describe the method used to per- 
form the relative calibration of the Pan- STARRS 1 data. 
First, in §2.1| we describe the general problem of the pho- 
tometric calibration of optical data. In §2.2| we describe 
the algorithm used to perform the relative cali bration of 
the survey, that of iPadmanabhan et ali ((2008). In §2.3| 
we then describe the details of the implementation of the 
algorithm for processing PS1 data. 

2.1. The Goal of Photometric Calibration 

The quantities of scientific interest in an imaging sur- 
vey are usually the astrometry and photometry of objects 
as a function of time. The photometry of an object gives 
the flux from that object reaching the earth within a 
filter bandpass. For linear detectors like the PS1 CCDs, 
instead the number of photons per second N reaching the 
detectors is directly measured, which is simply related to 
the flux / by 

N = Kf (1) 

if noise is ignored. The task of the photometric calibra- 
tion is to solve for the throughput K. 

Conventionally, the photometry of an object is given in 
magnitudes m, with m = — 2.5 lo g.f //o, where fn give s 
the AB magnitude reference flux (jOke and Gunnl fl983). 
Then Equation [T] becomes mi nst + Z = m, with the in- 
strumental magnitude given by mi nst = — 2.5 log TV and 
the zeropoint given by Z = 2.5 log K/q. We work in these 
logarithmic variables for the rest of this work, and so seek 
to determine the zero points Z of each observation in the 
survey. 

The zero points Z are determined by the light collect- 
ing efficiency of all of the components of the system. At 
a particular wavelength A, Z can be decomposed into a 
number of factors describing the system: 

Z = 2.5 log AT a T o T f T d f (2) 

Here A gives the collecting area of the telescope, and T a , 
T Q ,Tf, and Td the throughputs of the atmosphere, optics, 
filter, and detector, respectively. In principle, Z can be 
different for each star in each exposure, if, for instance, 
the filter throughput or detector efficiency varies over the 
focal plane. 

The relative calibration of the survey is concerned only 
with how Z varies from object to object in the survey. We 
can then separate Z into two terms, Z = Z a + Z r , where 
Z a is a constant giving the absolute zero point of the 
survey in a particular situation, and Z r gives the relative 
change in zero point from Z a . Unchanging components 
of Z, like the reference flux fo and mirror area A, affect 
only Z a and are independent of Z r . The model for Z r 
must be simple enough to allow its components to be 
constrained, but flexible enough to capture the variation 
in Z r . 

For a system like PS1, where the optical system, filter, 
and detector are essentially unchanged over the course 
of the night, we can suppose that A, T , Xy, and Td are 



constant over the course of the night. We can then encap- 
sulate the effect of all these terms on the zero point as a 
single term, a, for that night, with a = 2.5 log AT Q TfTd, 
and seek to measure how a varies over the course of the 
survey. 

We must also account for T a , the variation in the trans- 
parency of the atmosphere. We model the atmospheric 
transparency simply as 

2.5\ogT a = -kx (3) 

where k describes the effectiveness of the atmosphere at 
extin guishing light, and x 'is th e air mass of the observa- 
tion (Pa dmanabhan et q/.ll2QQ8[ ). The survey is executed 
so that all observations have low air mass; the largest air- 
mass included is 2.7 and the vast majority of images have 
air mass less than 1.6. Equation [3] strictly holds only for 
monochromatic light, or when the atmospheric extinc- 
tion does not vary with wavelength. This assumption is 
violated in the ypi band, where there are strong water 
absor ption features in t he atmospheric extinction (see $4] 
and lTonrv et ali (|2Q12l )). Moreover, k will only be con- 
stant over a night if the atmosphere is isotropic and un- 
changing in time, an assumption clearly violated when 
clouds are present. Still, we find that for most of the 
nights of the survey, the simpl e mo del 2.5 log T a = — kx 
is largely satisfactory (but see §4.21 for more details). 

We therefore present as a starting point for our photo- 
metric model the simple expression Z — a — kx, consis- 
tent with the above discussion. We ultimatel y ado pt a 
more complicated variation of this function in §2.31 The 
problem of photometric calibration then becomes to de- 
termine the parameters a and k of such a model for each 
night of the survey. We perform this calibratio n follow- 
ing the algorithm of IPadmanabhan et ali ([2008), finding 
the parameters that minimize the variance of repeat ob- 
servations of each star. 

2.2. Algorithm 

An optical survey provides instrumental magnitudes 
minst of objects in the sky. We may have several repeated 
observations of the same object, rai ns t,o,i with uncertain- 
ties cr 0> i, where o labels an object and i labels its obser- 
vations. We are ultimately interested in the calibrated 
magnitude m of the object, with m = mi nst + Z a + Z r: 
where Z a and Z r are defined in §2.11 We find Z r by 
minimizing 

^^(m ,,-m^) 2 /^ (4) 

o i 

where m is the average of m 0t i over all observations of 
object o. The absolute zero point Z a cancels out of this 
expression. 

Letting m be a vector with every observation of every 
object in the survey, and using a linear model for Z, 
then m = mi nst + Ap, where A is the design matrix 
for Z and p contains the parameters of the model for 
Z. For the simple model for Z described in £12.11 p is 
a vector containing a parameter a and k for each night 
of the survey. The vector m has length n Q b s , the total 
number of observations of all objects in the survey. The 
matrix A has dimensions 7i bs * ^pan where 72 par is the 
number of parameters in the model. Furthermore, let 
W be the n Q b s x ^obs matrix of weights, such that if mj 
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corresponds to an observation of object o, then 

^ Wijrrii = m (5) 

i 

We then want to solve 

= m - Wm (6) 

in a least-squares sense. Expanding m in terms of m inst , 
and rearranging, we obtain 

(1 - W)Ap = (W - l)m inst (7) 

We solve this in a least-squares sense using a simple diag- 
onal covariance matrix C, with diagonal elements equal 
to the photometric variance in measurement i. We im- 
pose an error floor of 0.01 mags on the photometric vari- 
ances, to prevent a few bright stars from dominating the 
fit. Letting A' = (1 - W)A and b = (W - l)m inst , this 
is an ordinary linear least squares problem with solution 

p = (A'TC^A'^A'C^b (8) 

We solve this equation directly to find p and hence Z r , 
the relative photometric calibration of the survey. 

The structure of A/ is i llustrated in more detail in 
iPadm anabh an et all (|2008f ). We mention one appeal- 
ing feature of the structure of the problem, however: 
if A', W, and mi nst are written as sums of A^, W$, 
and minst, i, where the terms in these sums contain only 
rows corresponding to observations of object i and are 
otherwise 0, then likewise A /T C _1 A / and A / C _1 b split 
into sums of terms involving only observations of a sin- 
gle object each. This simplifies the computation of the 
matrices A /T C -1 A' and A / C _1 b and allows the terms 
contributing to them to be computed in parallel without 
ever requiring the matrix A' to be built. This is critically 
important because A' has size n Q bs x w par , which for a 
survey of a billion observations and a model containing 
thousands of parameters contains trillions of elements. 
Because of the intrinsic parallelism of the problem, only 
A /T C _1 A' needs to be built, which is of manageable size. 
Likewise, the parallelism means that only the observa- 
tions of a single object, and not the observations of the 
entire survey, need to be simultaneously read into mem- 
ory. 

The matrix A /T C -1 A', computed as described above, 
is necessarily singular. The singularity occurs because 
inevitably, for any solution Z r , the solution Z r + c is 
equally good for any c. The constant c is degenerate with 
the absolute calibration of the survey. Depending on the 
model for Z r , other singularities may exist. Accordingly, 
we perform a singular value decomposition of A /T C _1 A'. 
Eigenvectors of A /T C _1 A' below a c ertain threshhold are 
fixed using priors, as described in Padmanabhan et al 
(2008). The singularities affecting the PS1 calibration 
are described in §2.31 

This method correctly derives best-fit parameters p de- 
scribing the relative photometric calibration. The above 
discussion has described the photometric calibration al- 
gorithm in general terms. We now describe the details 
of the calibration as applied to the PS1 survey. 

2.3. Details for Pan-STARRSl 

The most important question to be answered in apply- 
ing the photometric calibration algorithm to a particular 



TABLE 2 
Flat Field Seasons 



Season 


Begin Date 


End Date 


Season 1 




10 April 2010 


Season 2 


10 April 2010 


11 May 2010 


Season 3 


11 May 2010 


11 April 2011 


Season 4 


11 April 2011 





Note. — The dates marking the bound- 
aries of the different time periods for which 
independent flat fields are used in the 
photo metric model, from iFinkbeine r et al\ 
( 2012). Season 1 includes all data taken be- 
fore 10 April 2010, while season 4 includes 
all data taken after 10 April 2010. 

survey is the choice of model for Z r . The simplest rea- 
sonable model uses Z r — Ojti ky^x^ where a n describes 
the throughput of the optics, filters, and detector on a 
night n, and k n describes the transparency of the atmo- 
sphere on that night, as described in §2.11 We need to 
refine this model slightly for PS1. 

In wide-field surveys, frequently the system response 
depends on the position of the star in the focal plane, for 
exam ple, because of th e need for an illumination correc- 
tion (jHogg et a/.ll2Q01f ). The raw Pan-STARRSl images 
are already corrected for non-uniform detector through- 
put across the field of view according to a static flat 
field derived from the combination of dome flats and 
stellar photometry. We nevertheless solve for a new 
flat field using the wealth of data taken since the be- 
ginning of the survey. Comparison between SDSS and 
PS1 data indicated that the fla t field changed abruptl y 
three times during the survey (jFinkbeiner et all 12012). 
We accordingly fit for separate flat fields fij for the four 
different "seasons" to account for this behavior, so that 
Z = a n — k n x + fij. Here i indexes over locations in 
the focal plane, which we take to be the four quadrants 
of each of the 60 CCDs in the PS1 focal plane, and j 
indexes over the four seasons (Table EJ). 

Analysis of Z derived for images taken as part of the 
medium-deep survey revealed that the amount of flux 
registered by the PS1 system systematically varied de- 
pending on the image quality of the individual images, 
in the sense that flux is lost for images with a very small 
or, especially, very large FWHM. To account for this 
variation, we include in our model for Z a quadratic w 
in FWHM, leaving us with the model 

Z = a n - k n x + fa + w(F) (9) 

where F is the FWHM of the image in which the observa- 
tion was made. This model is summarized in Table [3j and 
is the final model we adopt to derive the Pan-STARRSl 
zero points. 

We have experimented with fitting only a single k term 
for the entire survey, rather than fitting one k term for 
each night. However, the best fit values of k can vary 
from night to night, by about 0.05 mags/airmass ( §4.3j) . 
That said, the median standard deviation in airmass x 
for all observations on a given night is only about 0.1, 
leading to an induced uncertainty of about 5 mmag from 
ignoring the variation in k. Nevertheless, because the 
extreme edges in declination of the survey, the north ce- 
lestial pole and declination —30°, must be observed at 
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TABLE 3 

Parameters of the Photometric Model 



Parameter 


Number 


Note 


a 


~ 200 


system (nightly) 


k 


~ 200 


atmosphere (nightly) 


f 


4 x 60 x 4 


illumination correction 


w 


2 


FWHM correction (quadratic) 



Note. — The parameters of the photometric model used in 
this work. The calibration is performed independently in each of 
the 5 Pan-STARRSl filters. Observations have been performed on 
about 200 nights in each filter, though the exact number ranges 
from 293 in zp\ to 190 in ypi. The illumination correction has 
one parameter describing each of the four quadrants of the 60 PS1 
CCDs, over four time periods. The constant term of the quadratic 
in w is not fit, because it is completely degenerate with a. 

relatively high air mass, neglecting the variation in k leads 
to errors in the photometric calibration of the survey on 
large angular scales. For this reason, we fit a k term for 
each night of the survey. 

Other models for Z could be adopted. In principle 
we have enough information to fit a zero point for every 
Pan-STARRSl image independently, though this would 
greatly diminish the stability of the solution on large an- 
gular scales. However, the Pan-STARRSl system has 
proven to be remarkably stable over the course of a night 
( §4.1)) . removing the need for a more finely grained cal- 
ibration, and leading us to adopt a simple model fitting 
only two parameters per night. 

For the photometric calibration we use only observa- 
tions taken at times when we believe the night to be pho- 
tometric. We define "photometric" here to mean that the 
simple model for Z described above is satisfied on that 
night to within about 20 mmag. In cloudy weather the 
Z behave erratically; when we find evidence for clouds 
greater than 20 mmag in Z, we manually flag that por- 
tion of the night as non-photometric. We flag about 25% 
of the images taken as non-photometric. We are able 
to recognize exposures as taken in non-photometric con- 
ditions only when they are discrepant with other over- 
lapping Pan-STARRSl exposures, or when they overlap 
the SDSS. As the Pan-STARRSl survey continues, our 
ability to flag non-photometric exposures will improve, 
though the analysis of §3] indicates that our performance 
is already good. 

We wish to use only secure observations of typical stars 
in the calibration, to avoid any bias in the calibration 
from anomalous measurements. We therefore use only 
objects for which at least one detection had estimated 
uncertainty less than 30 mmag. We also exclude any 
objects for which any measurement of that object had 
an instrumental magnitude less than —14.25, to avoid 
any detector saturation effects. We finally exclude any 
detections on images with FWHM greater than four arc- 
seconds. 

We use techniques for measuring mi nst that assume 
that the object being measured is a point source (i.e., 
PSF mags). Accordingly, we wish to include only point 
sources in the calibration, excluding the galaxies. For 
this purpose, we exclude any object for which more than 
25% of the detections of that object have PSF magnitude 
minus aperture magnitude greater than 0.1 mags. We 
also exclude any objects for which more than 10% of the 
detections of that object have m > c, with c equal to 19 



mag in gpi and rpi, 18.75 in zpi, and 18 in zpi and ypi. 
These correspond approximately to the magnitudes at 
which the SDSS finds that the number density of galaxies 
exceeds the number density of stars, and so these cuts 
further reduce the galaxy contamination in our analysis. 
Our tests indicate that resulting selection of stars is very 
clean, and varying the star-galaxy separation negligibly 
affects our results. 

We iterate the photometric calibration several times 
and clip discrepant observations and images on each iter- 
ation, reducing the clipping threshhold by a factor of two 
on each iteration until reaching 3<r, as described below. 
This reduces the sensitivity of the algorithm to outliers. 
In clipping detections, we compute A = m — fn for each 
detection. For each image z, we find the mean ji^ and 
standard deviation a a of A on that image. An image is 
clipped if /iA or a a is inconsistent with their respective 
distributions for all images at the clipping threshhold. 
About 3% of images are clipped in this process; most 
periods of non-photometric conditions have already been 
flagged by hand. A detection on an image is clipped if 
|A/cr| is greater than the clipping threshhold for that 
detection, where a is the sum in quadrature of the pho- 
tometric uncertainty for that detection and <ja for that 
image. About 4% of detections are clipped. These in- 
clude detections with problematic photometry and vari- 
able stars. 

We modified the algorithm of Padmanab han et all 
(2008) slightly to incorporate including an external 
source of photometry in the calibration. When enabled, 
if external photometry for an object is available, we set m 
for that object to be that given by the external photom- 
etry, so that for these objects the best fit solution simply 
minimizes the difference between the external and inter- 
nal photometry. We can then add, for instance, all of 
the photometry from the SDSS matching PS1 objects 
into the photometric calibration. We have experimented 
with both including and excluding the SDSS photometry. 
In the final photometric calibration we currently use for 
PS1, we include the SDSS photometry, in order to im- 
prove our ability to detect non-photometric conditions. 
The derived zero points vary only slightly depending on 
whether or not the SDSS is included as an external ref- 
erence in the calibration, by about 5 mmag rms £13.21 

We impose Gaussian priors on the parameters of the 
photometric model when the photometry poorly con- 
strains them, as described in Table |U We take the 
set of poorly constrained eigenvectors of the solution 
and find the best fit parameters for those eigenvec- 
tors such that the p riors are satisfied, as described in 
Padmanabh an et all (|2008f ). The priors serve three pri- 
mary functions. The first is to set the absolute calibra- 
tion of the survey; the second is to resolve the degeneracy 
between the a and / terms; and the third is to resolve the 
degeneracy between the a and k terms on nights when 
the range of airmasses probed is small. 

We impose a prior for the absolute calibratio n of the 
survey based on the work of iTonry et all (j2012l ). which 
uses photometry from standard stars and the Hubble 
Space Telescope to determine the absolute calibration 
of Pan-STARRSl. We choose a prior for the a terms 
to agree with the results of that work. When using 
the SDSS as a reference, we im pose this pr i or by us- 
ing the color transformations of iTonrv et all (j2Q12h to 
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TABLE 4 
Priors 



Parameter 


Prior 


cr 


Parameter 


Prior 


cr 


a g 


24.408 


1 


kg 


0.147 


0.05 


a r 


24.679 


1 




0.085 


0.03 


a t 


24.556 


1 


ki 


0.044 


0.02 


a z 


24.218 


1 


k z 


0.033 


0.02 


a y 


23.237 


1 


ky 


0.073 


0.03 


f 





0.02 


W 





0.1 



Note. — The Gaussian priors for the parameters in the pho- 
tometric model. The priors on w are irrelevant as w is always well 
constrained by the data. The priors on / serve only to break the 
degeneracy between a and /. The priors on a set the absolute 
calibration of the survey, an d are tuned to match the SDSS with 
the color transformations of Tonrv et al. (2012). The priors on k 
constrain k on nights when the range of airmasses probed by the 
survey is small, causing k and a to be degenerate. 

transform the SDSS magnitudes into the Pan-STARRSl 
bandpasses, and then the absolute calibration is fixed di- 
rectly by the photometric solution by reference to the 
color-corrected SDSS. 

The a and / terms in the photometric model have an 
exact degeneracy, in that Z is unchanged if all of the 
a are increased and / is decreased. This degeneracy is 
removed by enforcing a prior on all of the / to be zero, 
with an uncertainty of 20 mmag. 

The a and k terms are degenerate if the range of air- 
mass probed on a night is small. This degeneracy is obvi- 
ously exact on a night consisting of a single image. Then 
kx is just a constant for the night, and so any change in 
a can be canceled with a change in k. These degeneracies 
are removed by setting a prior on k to be equal to the 
value of k on a typical night. 

After we obtain an iterated photometric solution, we 
perform a final adjustment to the derived Z for each im- 
age. We can robustly compute ha,i fc> r each image i in 
the survey. Ordinarily we use fiA,i to clip images that 
are discrepant from the rest of the photometric solution. 
Having obtained a photometric solution, however, we use 
the /iA,i to improve the calibration of the survey by ad- 
justing the photometric solution Zi for image i by /iA,i- 
This induces no changes in the photometric calibration 
on large spatial scales, but does clean up the light curves 
of bright objects slightly. The standard deviation of /iA,i 
is about 5 mmag (§3j). 

3. RESULTS 

We compute the PS1 photometric calibration using 
about one billion observations in each filter, solving for 
about a thousand parameters for each filter. These pa- 
rameters describe the system zeropoint a, atmospheric 
transparency /c, and flat field, as well as an image qual- 
ity correction as described by Equation [9] We perform 
the computation in parallel over the available cores in 
our system, using the inherent parallelization described 
in §2.21 The computation time is dominated by reading 
data from disk for processing, and takes about two days. 
The parallelization of the computation and the database 
operations more generally were g reatly simplified by th e 
Large Survey Database software (|Juricllin preparationh . 

We find that the results of the photometric calibra- 
tion account for variations in the mean zero point per 
exposure of the Pan-STARRSl system to better than 10 
mmag rms, albeit with some areas of worse calibration. 



We demonstrate this accuracy, checking the internal con- 
sistency of the solution ( §3.1|) . the consistency with the 
SDSS ( §3.2p . and the consistency of the colors of stars in 
different parts of the sky ( §3.3p . 

3.1. Internal Consistency 

We first test the internal consistency of the photomet- 
ric calibration in three ways: by using simulated data, by 
examining the residuals from our photometric model over 
the sky, and by comparing with an alternative method 
for calibrating the MD fields. 

Tests with simulated data indicate that errors due to 
the statistical uncertainty of the observations are negli- 
gible. In these simulations, which do not include clouds, 
we take the actual set of observations of all of the stars 
used in this analysis. For each star, we then declare 
its true magnitude for the purposes of the simulation to 
be equal to the mean of all measurements of that star. 
We then convert these true magnitudes to instrumental 
magnitudes, using our adopted model for the zero points, 
with parameters drawn at random from the distributions 
given in Table HI Noise is then added to these simu- 
lated measurements consistent with estimates from the 
actual Pan-STARRSl data, to produce a set of simulated 
Pan-STARRSl observations. 

These simulated observations are then photometrically 
calibrated, and the recovered zero points for each ob- 
servation are compared with the true zero points. The 
standard deviation in the difference between the recov- 
ered magnitudes and true magnitudes is less than 0.2 
mmag when using a simulated external reference catalog 
designed to match the SDSS. If this simulated external 
reference catalog is not used, the standard deviation re- 
mains less than 0.2 mmag, though a small number of 
observations have zero points off by as much as 6 mmag 
in the ypi band. 

These simulations do not include systematic errors 
caused by deviations in the system throughput from our 
model (for instance, clouds). However, they verify that 
the data have the signal-to-noise necessary to constrain 
the model, and that the algorithm accurately recovers the 
photometric model parameters when given good data. 

We perform an internal test of the size of the devia- 
tions from our photometric model using A = m — fa. 
For each exposure, we compute [1a and <ta, the mean 
and standard deviation of A for all observations on a 
single exposure. When an individual exposure has a zero 
point inconsistent with our photometric model, ha will 
depart from zero. Exposures with highly variable point 
spread functions or other problems with the photometry 
will have large a a • The typical size of a a will depend on 
the brightness of the stars used in the calibration; here 
we are interested only in verifying that the a a are homo- 
geneous over the survey. Figure [2] shows a map of \/jLa\ 
and a a in each of the PS1 bands, using a PS1 calibration 
that does not include the SDSS as an external reference. 
The maps of \/jla\ are fairly uniform over the sky; a few 
isolated areas are slightly problematic, but no large scale 
trends are obvious and the Galactic plane makes no ap- 
pearance. The maps of a a are somewhat patchier and 
the Galactic center appears as an area of only slightly in- 
creased a A: a testament to the performance of the IPP in 
crowded fields. The a a have the tendency to be larger in 
areas with more observations, despite correction of a a by 
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y/N/ (N — 1), where N is the number of observations of 
the sky. This is likely because exposures taken as part of 
TTI pairs are especially photometrically consistent; the 
two exposures are taken only fifteen minutes apart with 
a common position angle, and so each detection lands on 
nearly the same pixels in the focal plane, and the point 
spread function has had little time to change. 

Figure [3] shows the same information as in Figure [2] 
as a histogram, giving the distribution of //a and a a- 
Unsurprisingly, JTX is near zero. The scatter in //a is 
only about 3 mmag. This is a lower bound on the actual 
uncertainty in our photometric calibration. The <ta are 
about 12 mmag, except in ypi, where the photometric 
scatter is about 17 mmag. 

The scatter in //a will be lower than the true uncer- 
tainty in the calibration because many stars have only 
been observed a few times, and often as part of corre- 
lated TTI pairs. We can eliminate this problem by lim- 
iting the exposures used to compute /iA to exposures on 
MD fields. In this case, the scatter in /iA goes up to 5 
mmag, while <ja remains about the same. This scatter 
in /iA is more realistic, and represents the accuracy of 
the photometric calibration that we could expect after 
covering the sky many times. However, sparsely covered 
portions of the sky may have photometric calibration er- 
rors much larger than found in these well-covered MD 
fields. 

As a final internal consistency check, we compare the 
zero points generated by the photometric calibration 
algorithm with t he ze ro points for the MD fields of 
iFinkbeiner et all (j2012[ ). In that work, a very loose pho- 
tometric model is adopted, which allows completely in- 
dependent zero points for each exposure and a separate 
flat field for each night. We find that the zero points of 
this work agree with those of IFinkbeiner et Til] (|2012l ) to 
about 5 mmag, consistent with our expectations from our 
internal tests of the model residuals on the MD fields. 

3.2. Consistency with the SDSS 

We can also verify the results of the photometric cal- 
ibration by comparing our zero points with those we 
would derive by forcing the photometry to match an ex- 
ternal reference as closely as possible. The SDSS has 
observed about one third of the sky, about half of the 
area that Pan-STARRSl has observed. We find zero 
points ^sdss for individual PS1 exposures that overlap 
the SDSS. We compare these zero points with the zero 
points obtained from the photometric calibration algo- 
rithm, performed without using the SDSS as a reference. 

We compute Zsdss by transforming the SDSS mag- 
nitudes of stars onto the Pa n-STARRSl s ystem using 
the color transformations of iTonry et all ()2012l ) . For 
each PS1 exposure, Z S bss = (^sdss - m in st). We use 
only SDSS stars with msDSS < 18 in the computation of 

^SDSS- 

The results of the comparison are shown in Figure HI 
The internal Pan-STARRSl zero points agree with the 
SDSS-based zero points to about 10 mmag in all bands, 
ranging from 7 mmag in rpi to 13 mmag in yp±. There 
are slight offsets in the mean zero point between the inter- 
nal and SDSS zero points. These means are determined 
by the absolute calibration of the survey, and are fixed 
by a prior; they provide no information about the rela- 
tive calibration that is the focus of this work. A small 



number of large photometric outliers do exist, and may 
be removed as the number of overlapping observations in 
the survey increases. 

The spatial structure of the differences between the in- 
ternal and SDSS-based zero points is of particular inter- 
est. Figure [5] shows maps of the mean difference between 
the calibrated PS1 magnitudes of objects used in the pho- 
tometric calibration and their color-transformed SDSS 
magnitudes, in pixels 0.2° on a side. The maps clearly 
show signs of errors in both the PS1 and SDSS photomet- 
ric calibration, suggesting that a simultaneous PS1-SDSS 
calibration would be valuable. The SDSS scan pattern 
is visible as narrow 3° stripes, approximately in right as- 
cension, while the PS scan pattern is seen as rectangles 
in right ascension and declination. The largest problems 
with the SDSS involve runs poorly connected to the main 
body of the SDSS; runs around a = 300° differ between 
PS1 and the SDSS by about 40 mmag. The internal PS1 
calibration shows clear ~ 25 mmag problems in gpi at 
(170°, 5°) and in zpi at (25°, 10°), to name a few, prob- 
ably due to unrecognized cloudy weather. Moreover, in 
zpi and ypi especially, parts of the maps are mottled at 
the 3° scale of the PS1 focal plane, indicating a problem 
with the PS1 photometry. Despite these problems, the 
overall agreement between the two sets of measurements 
is remarkable; the rms of these maps is about 11, 10, 11, 
12, and 16 mmag in grizypi. The results presented in 
Figure 3] are slightly better than these, because for that 
figure zero points were calculated for individual PS1 im- 
ages, averaging over multiple SDSS runs and over the 
PS1 photometric nonuniformity in the focal plane. 

3.3. Consistency of Stellar Colors 

We can additionally check the accuracy of the photo- 
metric calibration by testing the consistency of the stellar 
locus over the sky. This technique has the advantage over 
comparison with the SDSS that it can be applied over the 
entire sky. However, the presence of dust and intrinsic 
variations in stellar populations can cause the colors of 
the stellar locus to vary, limiting the effectiveness of this 
technique. 

We measure the color of main-sequence turn-off stars 
over the sky to test the consist ency of the st el lar lo - 
cus, following the technique of Schl afly et all j2010). 
This is one of a number of related techniques; we could 
alternatively have used the principal color analysis of 
Ivezic et ~au (120041 ) or the stellar locus regression of 
High et al\ (|2009[ ). Figure [6] shows maps of the color of 
the blue tip of the stellar locus in the PS1 bands grizypi. 

The dominant signal in Figure [6] is clearly from the 
inters tellar dust. After rem oving the dust accord- 
i ng to ISchlegel et all (|1998[ ) and lSchlaflv and Finkbeinerl 
(2011), at high Galactic latitudes most of the signal 
comes from problems with the photometric calibration. 
The most egregious example is in the ypi band at (25°, 
-25°), of 40 mmag. Over most of the sky, calibration 
errors are consistent with expectations from the compar- 
ison with the SDSS, ~ 15 mmag. 

4. DISCUSSION 

We interpret these results in the context of the stability 
of the Pan-STARRSl system. The system stability di- 
vides naturally into four different com pone nts: the pho- 
tometric stability over a single night ( §4.1|) . the stability 
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Fig. 2. — Maps of |/xa| (left panels) and cta (right panels) in grizypi. The x-axes give right ascension and the y-axes give declination, 
both in degrees. These give the consistency of the zero points and the scatter in the photometry over the sky. White to black is 0-10 mmag 
for the left panels and 5-25 mmag for the right panels. The Galactic center is barely visible as a region of increased cta near (266°, —29°). 
The <ta are substantially larger in the ypi band than in the other bands. 



o f the detector and optics over the course of the survey 
( 34.2|h the stability of the atmosphere over the course of 
the survey ( §4.3|) . and the stability of the flat field over 
the survey 

4.1. Nightly Photometric Stability of PS1 

The photometric model we have adopted to calibrate 
the survey assumes that the throughput of the system 
and atmosphere do not vary substantially over the course 
of a night. This assumption is occasionally violated, lead- 
ing us to remove about 25% of the nights on which sur- 
vey data are taken. The stability of the system on the 



remaining nights is excellent. Figure [3] showed that the 
typical model residuals are less than 5 mmag. However, 
the nightly stability can be shown more explicity in a plot 
of model residual as a function of time of night. Figure 
[71 shows a variety of plots of data taken on February 13, 
2010, during the first month of full science operation of 
the survey. The first panel shows a density plot giv- 
ing the distribution of A in mmag for each star in each 
image as a function of the hour in the night, with con- 
tours marking the mean and ±lcr of A for each image. 
Dashed blue lines mark ±20 mmag. The crosses indicate 
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Fig. 3. — Histograms of ha (left) and <xa (right) in grizy^i- The mean h and standard deviatio n a o f /xa and <ta are labeled for each 
filter. The scatter in represents a lower bound on the uncertainty in our photometric calibration f 33.1) 1. The a a are quite similar among 
the various bands, except in ypi, which has 50% more scatter than the other bands. 



the fraction of the detections clipped from each image in 
the photometric calibration, and their color indicates the 
band that each image was taken in, using the same colors 
as in Figure [8j The photometric stability on this night is 
very good; the majority of images have \ii/s\ < 5 mmag. 
The gpi images taken around 14 hours, however, are dis- 
crepant at the 20 mmag level; this is the source of the 
clear Pan-STARRSl calibration error evident at (170°, 
10°) in the g band in Figure [5j and seems to be one of 
the largest calibration errors remaining in the data. 

4.2. System Stability 

The scatter in the zero-point residuals from the simple 
nightly photometric model we adopt is about 5 mmag, 
showing that the system is stable over the course of a 
night. However, this stability is not simply a nightly 
phenomenon; the system has been very stable over the 
whole course of the survey. 

Ideally, we would measure the stability of the system as 
the scatter in the a-terms derived by the photometric cal- 
ibration. This procedure can overestimate the true scat- 
ter in the system throughput, however, because of the 
degeneracy between the a- and /c-terms on nights when 
the range of airmass probed by the survey is small. The 
median standard deviation of airmass in images taken in 
a single filter on a single night is only 0.1, rendering a and 
k substantially degenerate on most nights. Accordingly, 
we instead test the stability of the system by looking at 
the scatter in the zero points Z at 1.2 airmasses, which 
is close to the average airmass of the survey. We ignore 
the dependence of Z on w and /; the former could con- 
tribute slightly to the scatter in Z, but should play only 
a minor role, while the latter is required to have mean 
zero. This combination of a and k is well-constrained, 
but includes variation both from a and from kx. The 
scatter from kx at 1.2 airmasses is Xfla^. We find in §4.31 
that (Tk is less than about 0.05, so the atmosphere should 
contribute less than 60 mmag scatter to the zero points. 



Figure [8] shows the zero points a — kx at airmass 
1.2 for each photometric night of the survey in the 5 
Pan-STARRSl bands. The yp\ bands have been offset 
by 0.7 mags for clarity. The zero points are extraor- 
dinarily stable, except in yp\. In the other bands, the 
intrinsic throughput of the Pan-STARRSl system has 
varied by less the 20 mmag rms. A change in the system 
throughput near Modified Julian Day 55524 of 20 mmag 
is obvious in griz. This stability is surprising given that 
we expect as much as 60 mmag scatter in zero point from 
the airmass term alone, and requires that the scatter in 
k is actually less than 0.02 in grizp\. 

The ypi band has by far greater scatter in its nightly 
zero points than the other bands. This is presumably 
the result of two factors, both owing to the presence of 
strong water vapo r absorption bands that overlap ypi 
(jTonrv et aZ.I l2Q12h . First, because the atmospheric ab- 
sorption is not approximately constant over ypi, model- 
ing the ypi zero point simply as linear in airmass is not 
appropriate. Second, the depth of the absorption bands 
varies with the amount of precipitable water vapor in the 
air, leading to varying zero points. An improved pho- 
tometric calibration would then include more freedom 
in the ypi zero points, potentially incorporating addi- 
tional information about water vapo r in the air, as has 
been proposed in a number of works dStubbs et all 120071 : 
iBurke et all 120101 : iBlake and Shawl |2011|), and adding 
color-airmass and color-water vapor terms to the pho- 
tometric model. 

The work of iTonrv et all (|2012[ ) considers the depen- 
dence of zero point on airmass in more detail. In partic- 
ular, that work finds that adopting 



2.5 log T a = -kx 



0.68 



(10) 



more accurately describes the relation between zero point 
and airmass in the ypi band. However, we find that 
owing to the small range of airmass probed on a single 
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Fig. 4. — Comparisons between the zero points of this work and zero points derived relative to the SDSS, for the niters grizypi (rows). 
The left column gives the zero points of this work relative to the zero points derived from the SDSS, both in magnitudes. It also gives 
the number N of images used in the comparison. The middle column gives the difference in mmag between the two zero points with the 
SDSS zero points, in magnitudes. The third column gives a histogram of the differences in mmag, along with their mean fi and standard 
deviation a. Unrecognized non-photometric conditions lead to points with low SDSS-derived zero points and large positive differences 
between SDSS- and PS-derived zero points. 



night, adopting this relation alters the zero points de- 
rived in this work only negligibly (~ 1 mmag). Even so, 
we note that extrapolating the PS1 magnitudes of stars 
to the magnitudes that would be observed above the at- 
mosphere requires taking the nonlinearity into account. 
That problem is one of absolute calibration; in this work 
we consider only the relative calibration of the survey. 

4.3. Atmospheric Stability 



The stability of the atmosphere over Haleakala places a 
fundamental limit on the stability of the Pan- STARRS 1 
system. The atmospheric transparency, however, can 
only be separated reliably from the system throughput 
when a wide range of airmass is probed over the course 
of a night. On the typical night, however, the standard 
deviation in airmass of all observations in a single filter 
is less than 0.1. In the zpi and zpi bands, this corre- 
sponds to less than 5 mmag of atmospheric extinction, 
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Fig. 5. — Maps of the difference between the color-corrected SDSS magnitudes of stars and the internally-calibrated Pan-STARRSl 
magnitudes of the same stars in the filters grizypi (rows). The x-axes give right ascension and the y-axes give declination. The rms of 
the maps is about 10 mmag. Narrow stripes in right ascension are symptomatic of problems with the SDSS photometric calibration, while 
rectangles in right ascension and declination indicate problems with the PS1 calibration. The filter used for each map is indicated in the 
lower left, while the rms of the map is indicated in the lower right. 
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Fig. 6. — Maps of the color of the main-sequence turn-off (MSTO) stars observed by Pan-STARRSl. The four rows give the colors g — r, 
r — i, i — z, and z — y. T he l eft column gives the obse rved MSTO color, while the right column gives the color corrected for dust according 
to Schlegel et al. (1998) and Schlaflv and Finkbeiner (2011). Clearly the dust accounts for most of the signal in these maps, but problems 
with the photometric calibration are also evident, for example, in y-pi at (25°, —25°). Blank areas have no Pan-STARRSl observations in 
photometric weather in one of the two bands making up the relevant color. Black and white are ±0.1 mag of the median color of each map. 



comparable in size to the typical model residuals, mak- 
ing it impossible to reliably measure the airmass term k 
in the photometric model for that night. 

The problem is further complicated by the fact that the 
photometric calibration has access only to the differences 
in magnitudes of observations of the same stars. Even 
when a night covers a wide range of airmass, if on another 
night the same stars were observed at the same air masses, 
the calibration could reliably determine that those nights 



had the same fc-term but not what that fc-term was. 

We use two techniques to test the variation in k from 
night to night. We limit ourselves to nights where cr x , the 
standard deviation in airmass examined on that night, is 
greater than 0.1, and then look at the standard deviation 
in k for those nights. We determine k from our photo- 
metric model, as wel l as using only zero points derived 
from the SDSS ( Q . 

Measurements of fc-terms derived from PS1-SDSS com- 
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Fig. 7. — The performance of the Pan-STARRSl photometry on 13 February 2010. The figure gives the distribution of residuals A, in 
mmag, of magnitudes of stars as observed on this night from the mean magnitudes of these stars, as a function of hour during the night. 
Each column in the figure corresponds to an image. The contours give the mean and zblcr of A for each image. Crosses give the fraction 
of observations in each image clipped from the calibration, with 0% corresponding to the bottom of the plot and 100% the top of the plot. 
The colors of the crosses give the band the image was taken in; the colors are the same as in Figure [8] The mean of A for each image is 
small throughout the nig ht, e xcept for a 20 mmag deviation around hour 14. This is one of the largest remaining photometric calibration 
errors in the survey. See ^4 II for details. 

We can repeat this analysis using /c-terms derived di- 
rectly from the photometric calibration. Because the PS1 
survey region covers a larger range of declination than the 
SDSS covers, this allows more high-airmass observations 
to be included. Using nights with a x > 0.1 indicates that 
cr k is 0.02, 0.02, 0.02, 0.04, and 0.05 in grizy P1 . Using 
only nights with <j x > 0.2 gives a k equal to 0.02, 0.02, 
0.02, 0.02, and 0.04, which is consistent with our con- 
straint from £14.21 If a k in y P \ is in fact about 0.04, then 
both the Pan-STARRSl system and atmosphere are less 
stable in y P1 than in grizpi. 
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Fig. 8. — The system zero point a — kx at airmass 1.2 derived 
by the photometric calibration for each night of the Pan-STARRSl 
survey. The mean \i (in mag) and standard deviation a (in mmag) 
of the zero points is given for each band. The ypizero points have 
been offset by 0.7 mags for legibility. The intrinsic system stability 
is better than 20 mmag in griz. The ypi band is by far less stable 
than the other bands, with 61 mmag of scatter in zero point, pre- 
sumably owing to the sensitivity of ypi to water in the atmosphere 
(see JO]). 

parisons on nights with <j x > 0.1 indicate that a k is 
0.03, 0.04, 0.03, 0.07, and 0.07 in grizy P1 . These re- 
sults are sensitive to the limit on a x \ taking <j x > 0.2 
results in a k equal to 0.03, 0.03, 0.03, 0.04, 0.09, though 
only about ten nights of the survey have a x that large. 
Reducing <j x increases a k in all bands, presumably be- 
cause ct/c becomes dominated by slight deviations from 
the photometric model. These estimates of o~k all ex- 
ceed our estimates for a k based on the stability of the 
total throughput at 1.2 airmass ( §4.2j) . suggesting that 
the 20 mmag of scatter in zero points is dominated by 
atmospheric variations rather than variations in the PS1 
system, and that the SDSS-based a k are overestimates. 



4.4. Detector Stability 

The stability of the Pan-STARRSl detector can be 
tested independently from the atmosphere. The raw PS1 
images are flat-fielded with a single, static flat field de- 
rived from dome flats and stellar photometry taken early 
in the survey. The photometric calibration derives four 
independent flat fields for different time periods in the 
survey (Table [2]). The time variation in these flat fields 
tests the stability of the detector. 

Figure [9] shows the mean flat field and the difference 
from the mean flat field for the four flat field seasons 
in grizypi. Each image consists of 16 x 16 pixels, de- 
scribing the four quadrants of each of the 60 OTA CCDs 
composing the PS1 focal plane. The mean absolute dif- 
ference of the flat for each season from the mean flat 
field is less than 5 mmag in all of the filters. Likewise, 
the standard deviation of the differences of the flat fields 
from the mean is less than 6 mmag in all of the filters. 
The bright ring present in the flat fields from the first 
season indicates that the edge of the focal plane in the 
first season was too faint and required correction. This 
may be due to problems with the photometry in images 
with poor image quality, which was especially problem- 
atic early in the survey and around the edge of the focal 
plane. 

The mean flat has a scatter of about 10 mmag in all of 
the filters. Some striping in the mean flat suggests that 
one side of the PS1 OTA CCDs has a ~ 10 mmag dif- 
ferent throughput from the other side. These variations 
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Fig. 9. — The flat fields, in mmag, derived by the photometric 
calibration, in the grizypi filters. The first column gives the mean 
flat field for the survey in each filter, while the later columns give 
the difference between the flat field for that season (Table [2} and 
the mean flat field. The standard deviation in the mean flat is 
about 10 mmag, and the mean absolute residual of the difference 
between this mean flat and four flat seasons is less than 5 mmag, 
in all of the filters. 

should have been removed by the static PS1 flat field, 
yet are nevertheless robustly present in the photometric 
calibration and in comparison with the SDSS. This may 
point to an anomaly occuring during the observations of 
the stars used to construct the static PS1 flat field. 

5. CONCLUSION 

We present the photometric calibration of the first 1.5 
years of Pan-STARRSl survey data. The per-image zero 
points we measure agree with those computed indepen- 
dently relative to the SDSS with rms scatter of 8.1, 7.1, 
9.0, 11.5, and 12.7 mmag in grizypi. This indicates 
that both surveys have zero points accurate at that level, 
when averaged over the 3° field of view of PS1. On sev- 
eral arcminute scales, photometric nonuniformities over 
the PS1 field of view and striping in the SDSS start to 
contribute, but the rms scatter remains only 11, 10, 11, 
12, and 16 mmag. We anticipate that as the PS1 survey 
continues, the accuracy of this calibration will improve 
as repeat observations reveal slight deviations from the 
simple photometric model we adopt. Internal tests of 
the calibration indicate that we may be able to achieve 
overall image zero point accuracy as good as 5 mmag. 

This accuracy renders the photometric calibration a 
minor source of s ystematic error in the PS1 photometry. 
iFinkbeiner et all (|2012[ ) has discovered a non-linearity 
in the PS1 photometry that can bias the photometry of 
faint sources by a few hundredths when image quality is 
poor. Likewise, poor image quality in some PS1 images 
lead to systematics in the PSF magnitudes that dominate 
the photometric calibration errors. 

This photometric calibration demonstrates that 10 
mmag calibration accuracy is possible using the survey 
strategies ty pical of major up coming surveys, like the 
Sky Mapper (iKeller et aZl I2QQ7T ). the Dark Energy Sur- 
vey (Flausher 2005, DES), the Hyper Suprime-Cam Sur- 
vey (jTaka da 2010, HSC) , and the Large Synoptic Survey 
Telescope (jTvsonl I2002L LSST). These surveys operate 
in a mode much like PS1, repeatedly imaging the sky 
one filter at a time. This is in contrast to the SDSS, 



which operated in a drift-scanning mode and nearly- 
simultaneously imaged the sky in each of its filters. We 
note, however, that the simultaneous five-color imaging 
of the SDSS leads to the SDSS having better color accu- 
racy than magnitude accuracy, and that it makes detect- 
ing non-photometric weather easier. 

The calibration also demonstrates the possibility of 
photometrically calibrating the ypi band, which includes 
a strong water vapor absorption feature at about 940 nm. 
Despite the variability of this feature, we still achieve 
zero point accuracy of 13 mmag in this band. This fig- 
ure may be further improved by including the amount of 
precipitable water vapor in the air into the calibration. 
SkyMapper, DES, HSC, and the LSST all also intend to 
observe in filters including this same feature. 

The zero points achieved by this work show that the 
PS1 system is photometrically stable. Zero points ex- 
trapolated to 1.2 air masses every photometric night have 
rms scatter of less than 20 mmag in grizpi, over the 
course of the survey. Moreover, these 20 mmag of scatter 
are dominated by the variability of the atmosphere. The 
stability of the PS1 optical system and detector is par- 
ticularly impressive given the continuous improvements 
to the system over its first 1.5 years of operation, to re- 
duce sky brightness, improve image quality, and defeat 
camera artifacts. 

The level of calibration accuracy we have achieved will 
enable many PS1 survey goals. The discovery of satellites 
of the Milky Way, cosmological investigations of super- 
novae and galaxy clustering, and studies of interstellar 
reddening all require accurate photometric calibration. 
When the PS1 data become public, this calibration will 
provide another benefit to the community: a set of well- 
calibrated observations of stars covering most of the sky. 
Together with the absolute calibr ation of the PS1 data 
as described in lTonrv et all (|2012f h this work provides an 
absolutely calibrated set of standard stars over the entire 
sky north of declination —30°, going much fainter than 
current data sets. 
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